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Abstract 



o 

(N 

O The transverse structure of the nucleon as probed in hard exclusive processes plays 

critical role in the understanding of the structure of the underlying event in hard 
collisions at the LHC, and multiparton interactions. We summarize results of our 
recent studies of manifestation of transverse nucleon structure in the hard collisions 
t-H at the LHC, new generalized parton distributions involved in multiparton interactions, 

presence of parton fluctuations. The kinematic range where interaction of fast partons 
of the projectile with the target reach black disk regime (BDR) strength is estimated. 

O ! We demonstrate that in the BDR postselection effect leads to effective fractional energy 

' losses. This effect explains regularities of the single and double forward pion production 

in dAu collisions at RHIC and impacts on the forward physics in pp collisions at the 
LHC. 



1 Introduction 

^} The start of the LHC puts at a forefront the task of the describing the high energy 

00 proton - proton collision events in the whole their complexity. In particular, to search 

for new particles it is necessary to understand the structure of the underlying struc- 
ture of events with dijets. The knowledge of the inclusive cross sections of hard binary 

y—i collisions, which are expressed through the convolution of the hadron parton distribu- 

tion functions (PDF's) and the hard parton - parton interaction cross section, is not 
sufficient for these purposes. 
^ A natural framework for description of the complete picture of the high energy inter- 

K^j action is the impact parameter representation of the collision. Indeed, in high-energy 

pp scattering angular momentum conservation implies that the impact parameter b 
becomes a good quantum number, and it is natural to consider amplitudes and cross 
sections in the impact parameter representation. The nucleon's light -cone wave func- 
tions, describing the partonic structure at a low resolution scale, can be expressed 
in terms of the longitudinal momentum fractions of the partons and their transverse 
positions relative to the center-of-momentum: ip p (xi, pi). 
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Figure 1: pp collisions at small and large impact parameters. Transverse and side views. 
Dark (gray) disks correspond to the areas occupied by hard (soft) partons. 



The studies of exclusive hard processes lead to the conclusion that the transverse 
distribution of gluons with x > 10~ 3 than the distribution of soft partons involved 
in generic inelastic pp collisions (Sections 2 and 3). As a result, the events with a 
dijet trigger should occur, in average, at much smaller impact parameters than the 
minimum-bias inelastic events - see Fig. [T] Probability of multiple soft and hard 
interactions is much higher for head-on collisions than for peripheral collisions. Hence 
one expects a much more active final states for the dijet triggered events than for the 
minimum- bias events. 

To describe the transverse geometry of the pp collisions with production of a dijet 
it is convenient to consider probability to find a parton with given x and transverse 
distance p from the nucleon center of mass, fi(xi,pi). This quantity allows a formal 
operator definition, and it is referred to as the diagonal generalized parton distribu- 
tion(GPD). It is related to non-diagonal GPDs which enter in the description of the 
exclusive meson production (Sect. 2). The transverse geometry of the pp collision with 
production of a dijet is represented in Fig. 2. 

For inclusive cross section in the pQCD regime the transverse structure does not 
matter - the cross section is expressed through the convolution of parton densities. 
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Figure 2: Transverse geometry of hard collisions. 

Indeed, we can write 

<j h oc J d 2 bd 2 p 1 d 2 p 2 S(p 1 +b-p2)fi(xi,pi)f 2 (x2,p2)a'2->-2 = 

J (fbd 2 pid 2 p 2 fl(xi, pi)f 2 (x 2 , P2)&2^2 = fl(xi)fo(x2)(T2-+2- (1) 

Here at the last step we used the relation between diagonal GPD and PDF: J d 2 pfj (x, p,Q 2 ) = 
fAr.O-). 

At the same time, as soon as one wants to describe the structure of the final state in 
production of say heavy particles, it is important to know whether a dijet production 
occurs at different average impact parameters than the minimum-bias interactions. 
We will argue below that dijet trigger selects, in average, a factor of two smaller 
impact parameters. This implies that the multijet activity, energy flow should be much 
stronger in these events than in the minimum-bias events. Obviously, the magnitude 
of the enhancement does depend on the transverse distribution of partons and on 
correlation between the partons in the transverse plane. This information becomes 
available now. 

2 Transverse structure of the nucleon wave function 

The basis for the quantitative analysis of the transverse nucleon structure is provided by 
the QCD factorization theorem for exclusive vector meson (VM) production 1 which 
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states that in the leading twist approximation the differential cross section of the pro- 
cess 72 + p — > VM + p is given by the convolution of the hard block, meson wave 
function and generalized gluon parton distribution, g(x\,X2,t\ Q 2 ), where Xi,X2 are 
the longitudinal momentum fractions of the emitted and absorbed gluon (we discuss 
here only the case of small x which is of relevance for the LHC kinematics) . Of par- 
ticular interest is the generalized parton distribution (GPD) in the "diagonal" case, 
g(x, t\Q 2 ), where x\ = xi and denoted by x, and the momentum transfer to the nucleon 
is in the transverse direction, with t = — (we follow the notation of Refs. [HE]). 
This function reduces to the usual gluon density in the nucleon in the limit of zero 
momentum transfer, g{x,t = 0|Q 2 ) = g(x\Q 2 ). Its two-dimensional Fourier transform 

g{x lP \Q 2 ) ee J^e^g(x,t = -Al\Q 2 ) (2) 

describes the one-body density of gluons with given x in transverse space, with p ee \p\ 
measuring the distance from the transverse center-of-momentum of the nucleon, and is 
normalized such that / d 2 p g(x, p\Q 2 ) = g(x\Q 2 ). It is convenient to separate the in- 
formation on the total density of gluons from their spatial distribution and parametrize 
the GPD in the form 

g(x,t\Q 2 ) = g(x\Q 2 ) F g (x,t\Q 2 ), (3) 

where the latter function satisfies F g (x,t = 0\Q 2 ) = 1 and is known as the two-gluon 
form factor of the nucleon. Its Fourier transform describes the normalized spatial 
distribution of gluons with given x, 

F g (x, P \Q 2 )^ j^e i ^F g (x,t = -Al\Q 2 ), (4) 

with J d 2 p F g (x, p\Q 2 ) = 1 for any x. 

The QCD factorization theorem predicts that the t-dependence of the VM produc- 
tion should be a universal function of t for fixed x (up to small DGLAP evolution 
effects). Indeed the t-slope of the J/ip production is practically Q 2 independent, while 
the t-slope of the production light vector mesons approaches that of J/ip for large Q 2 . 
The independence of the measured differential cross sections of exclusive processes at 
\t\ < 1 GeV 2 is commonly described either by an exponential, or by a dipole form in- 
spired by analogy with the nucleon elastic form factors. Correspondingly, we consider 
here two parametrizations of the two-gluon form factor: 

f exp{B g t/2), 

F g {x,t\Q 2 ) = \ , , ^ 2 ( 5 ) 
i (l-t/m 2 g )- 2 , 

where the parameters B g and m g are functions of x and Q 2 . The two parametrizations 
give very similar results if the functions are matched at \t\ — 0.5 GeV 2 , where they are 
best constrained by present data (see Fig. 3 of Ref. [1]); this corresponds to [3] 

B g = 3.24/m 2 . (6) 
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The analysis of the HERA exclusive data leads to 

B g (x) = B g0 + 2a' g ln(x /x), (7) 

where x = 0.0012, B g0 = 4.1 (±g;|) GeV~ 2 ,a^ = 0.140 (±g;gf) GcV~ 2 for Q 2 , - 3 
GeV 2 . For fixed x, B(x,Q 2 ) slowly decreases with increase of Q 2 due to the DGLAP 
evolution j2]. The uncertainties in parentheses represent a rough estimate based on 
the range of values spanned by the HI and ZEUS fits, with statistical and systematic 
uncertainties added linearly. This estimate does not include possible contributions to 
a' g due to the contribution of the large size configurations in the vector mesons and 
changes in the evolution equation at —t comparable to the intrinsic scale. Correcting 
for these effects may lead to a reduction of a' g and hence to a slower increase of the 
area occupied by gluons with decrease of x. 

It is worth noting here that the popular Monte Carlo description of the pp collisions 
at the collider energies - PYTHIA uses x-independent transverse distribution of partons 
described by the sum of two exponentials. This distribution roughly equivalent to the 
dipole parametrization with m 2 « 2GeV 2 .6] which is hardly consistent with the data 
on J /ip photoproduction, see dashed line in Fig. 3. For smaller x the difference is even 
larger since the transverse size increases with decrease of x - see Eq. 7. 



3 Impact parameter distribution of proton— proton 
collisions with dijet production 

Using the information on the transverse spatial distribution of partons in the nucleon, 
one can infer the distribution of impact parameters in pp collisions with hard parton- 
parton processes [2] • It is given by the overlap of two parton wave function as depicted 
in Fig. [4] 

The probability distribution of pp impact parameters in events with a given hard 
process, P 2 (xi, X2, b\Q 2 ), is given by the ratio of the cross section at given b and the 
cross section integrated over b. As a result 

P2(x!,x 2 ,b\Q 2 ) = Jd 2 Pl Jd 2 p 2 S^(b-p 1+ p 2 ) 

xF g (x uPl \Q 2 ) F g (x 2lP2 \Q 2 ), (8) 
which obviously satisfies the normalization condition 

J d 2 bP 2 {x u x 2 MQ 2 ) = i- (9) 

This distribution represents an essential tool for phenomenological studies of the un- 
derlying event in pp collisions [H |3] • We note that the concept of impact parameter 
distribution is also used in MC generators of pp events with hard processes [3 H] , albeit 
without making the connection with GPDs, which allows one to import information 
on transverse nucleon structure obtained in the independent measurements. 
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Figure 3: t-dependence of the exclusive J/ip photoproduction data from the FNAL 
E401/458 experiment [5]. Solid line: t-dependence obtained with the exponential 
parametrization of the two-gluon form factor, Eq.5 (the slope of the J/ip cross section 
is Bj/^p = B g + AS, where AB « 0.3 GeV -2 accounts for the finite size of the J/ip] see [3] 
for details). Dashed line: i-dependence obtained with the dipole parametrization, Eq.5. 
Dotted line: t-dependence obtained with PYTHIA, effectively corresponding to a dipole 
form factor with m 2 ~ 2 GeV 2 . 

For the two parametrizations of Eq. Eq. Q leads to ( for x = x\ — X2) 

P 2 {x,b\Q 2 ) = { , , (10) 

{ [m 2 g l(l2v)](m a b/2y K 3 (m g b), 

where the parameters B g and m g are taken at the appropriate values of x and Q 2 . 

The derived distribution should be compared to the distribution of the minimum- 
bias inelastic collisions which could be expressed through T(s, b) that is the profile 
function of the pp elastic amplitude (r(s, b) = 1 if the interaction is completely absorp- 
tive at given b) 

p- m ( s ,b) = [i-|i-r( s ,6)| 2 ] Ax in ( s ), (ii) 

where J d 2 bP- m (s,b) = 1. 

Our numerical studies indicate that the impact parameter distributions with the 
jet trigger (Eq |10[ ) are much narrower than that in minimum-bias inelastic events at 
the same energy (EqpTj) and that 6-distribution for events with a dijet trigger is a 
very weak function of the px of the jets or their rapidities, see Fig. 5. For example for 
the case of the pp collisions at ^Js — 7 GeV we find the median value of b, b me dian ~ 
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Hard interaction 




Figure 4: Overlap integral of the transverse spatial parton distributions, defining the impact 
parameter distribution of pp collisions with a hard parton-parton process, Eq. Q. 

1.18 fm and b me dian ~ 0.65 fm for minimum-bias and dijet trigger events (Fig. 5b)[3]. 
Since the large impact parameters give the dominant contribution to <7meZ our analysis 
indicates that there are two pretty distinctive classes of pp collisions - large b collisions 
which are predominantly soft and and central collisions with strongly enhanced rate 
of hard collisions. We refer to this pattern as the two transverse scale picture of pp 
collisions at collider energies [3J. 

A word of caution is necessary here. The transverse distance b for dijet events is 
defined as the distance between the transverse centers of mass of two nucleons. It may 
not coincide with b defined for soft interactions where soft partons play an important 
role. For example, if we consider dijet production due to the interaction of two partons 
with x ~ 1, Pi,p2 *~ since the transverse center of mass coincides with transverse 
position of the leading quark in the x — > 1 limit. As a result b for the hard collision will 
be close to zero. On the other hand the rest of the partons may interact in this case 
at different transverse coordinates. As a result, such configurations may contribute to 
the inelastic pp cross section at much larger b for the soft interactions. However for 
the parton collisions at Xi,X2 <C 1 the recoil effects are small and so two values of b 
should be close. 

The present pp LHC data already provide important tests of this picture. Let 
us consider production of the hadron (minijet) with momentum px- The observable 
of interest here is the transverse multiplicity, defined as the multiplicity of particles 
with transverse momenta in a certain angular region perpendicular to the transverse 
momentum of the trigger particle or jet (the standard choice is the interval 60° < 
| A0| < 120° relative to the jet axis; see Ref. [3] for an illustration and discussion of the 
experimental definition). In the central collisions one expects a much larger transverse 
multiplicity due to the presence of multiple hard and soft interactions. At the same time 
the enhancement should be a weak function of px in the region where main contribution 
is given by the hard mechanism [5J [3] . The predicted increase and eventual flattening 
of the transverse multiplicity agrees well with the pattern observed in the existing 
data. At s/s = 0.9 TeV the transition occurs approximately at pr.crit ~ 4GeV, at 
^/s = 1.8 TeV at p TiCrit « 5GeV, and at p T , crlt = 6 - 8GcV for 7 TeV [HUE]- We 
thus conclude that the minimum px at which particle production due to hard collisions 
starts to dominate significantly increases with the collision energy. This effect is likely 
to be related the onset of the high gluon density regime in the central pp interactions 
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Figure 5: (a) Impact parameter distributions of inelastic pp collisions at yfs = 7TeV. 
Solid (dashed) line: Distribution of events with a dijet trigger at zero rapidity, 3/1,2 = 0, 
for pt = 100 (10) GeV cf. Eq. (10). Dotted line: Distribution of minimum-bias inelastic 
events,cf. Eq. ( 11 ). (b) Dependence of median b on pt for different rapidities of the dijets. 



since with an increase of incident energy partons in the central pp collisions propagate 
through stronger and stronger gluon fields. 

Many further tests of the discussed picture which are suggested in Ref. |3, will be 
feasible in a near future. They include (i) Check that the transverse multiplicity does 
not depend on rapidities of the jets, (ii) Study of the multiplicity at y < for events 
with jets at yi ~ J/2 ~ 2. This would allow to check that the transverse multiplicity 
is universal and that multiplicity in the away and the towards regions is similar to 
the transverse multiplicity for y < 0. (iii) Studying whether transverse multiplicity is 
the same for quark and gluon induced jets. Since the gluon radiation for production 
of W ± ,Z is smaller than for the gluon dijets, a subtraction of the radiation effect 
mentioned below is very important for such a comparisons. 

Note that the contribution of the jet the fragmentation to the transverse cone as 
defined in the experimental analyses is small but not negligible especially at smaller 
energies (yfs — 0.9TeV). It would be desirable to use a more narrow transverse cone, 
or subtract the contribution of the jets fragmentation. Indeed, the color flow contri- 
bution [12] leads to a small residual increase of the transverse multiplicity with px- 
However the jet fragmentation effect depends on pt rather than on yfs. Hence it does 
not contribute to the growth of the transverse multiplicity, which is a factor of ~ 2 
between yfs = 0.9TeV and yfs = 7.0 TeV. In fact, a subtraction of the jet fragmenta- 
tion contribution would somewhat increase the rate of the increase of the transverse 
multiplicity in the discussed energy interval. This allows to obtain the lower limit for 
the rate of the increase of the multiplicity in the central (< b >^ 0.6 fm) pp collisions 
of s . It is a bit faster than the s dependence of multiplicity in the central heavy ion 
collisions. 
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Figure 6: Geometry of two hard collisions in impact parameter picture. 



4 Multiparton interactions 

Probability of the multiparton interactions which result in two hard collisions grows 
rapidly with the increase of the incident energy. Understanding of such processes is 
important for detailed understanding of the high energy QCD dynamics as well as for 
practical purposes - estimating backgrounds for the searches for new particles. By 
exploring the difference in scales between soft and hard QCD processes and space-time 
structure of Feynman diagrams we derive within pQCD the general formulae for the 
two dijet production in pp collisions [13j [14] and find that it contains two contribu- 
tions. The contribution which dominates in a wide range of Xi is the 4 — > 4 process 
which matches the intuitive geometric picture depicted in Fig. [6] However, a consistent 
pQCD treatment requires that one also takes into account 3 — !• 4 double hard inter- 
action processes that occur as an interplay between large- and short-distance parton 
correlations [14]. Such contributions are not taken into consideration by approaches in- 
spired by the parton model picture. This contribution takes into account correlations 
between the partons induced by the pQCD evolution. Our analysis indicates that this 
contribution becomes important only for x < 10~ 3 . 

The 4 — > 4 cross section for the collisions of hadrons "a" and "b" has the form 1 13): 

d(Ti = J J^v2 J dxi J dx 2 J dx 3 J dx 4 



_>> da 13 da 2A 

x D a (x 1 ,x 2 ,p 1 ,p 2 , A)D b (x 3 ,x 4 ,pl,pl,-A)— li -^dt x di 2 - (12) 

dt\ dt 2 



Here D a {x\, x 2l p\ 1 p 2l A) are the new "double" GPDs for hadrons "a" and "b" 

D( Xl ,x 2 y 1 y 2 ,t) = f2 J d ^ 2 ^ 2 eiA-kl)e{ P l-kl) 



, 2 , „i i=n i=n 

IT / IT dx t (2n) 3 s(j2^~mJ2^ 

ij£l,2 ^ n ' •'° i^l,2 i=l i=l 
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_i_ — y — y — }■ — y -* 
xip n (x 1 ,k 1 ,X2,k 2} . } ki,Xi..)ipZ{xi,ki + ^,x 2 ,k 2 - A , x 3 , fc 3 , ...)■ (13) 

Note that this distribution is diagonal in the space of all partons except the two partons 
involved in the collision. Here ip is the parton wave function normalized to one in the 
usual way. An appropriate summation over color and Lorentz indices is implied. 
Within the parton model approximation the cross section has the form: 

(J 4 = (j 1 a 2 /TTRf nt , (14) 

where a± and a 2 are the cross sections of two independent hard binary parton inter- 
actions. The factor t characterizes the transverse area occupied by the partons 
participating in the two hard collisions. It also includes effect of possible longitudinal 
correlations between the partons. 



Eq. (12 1 leads to the general model independent expression for 



——^D(xi,x 2 , —~A)D{x 1 ,x 2 ,li), (15) 

in terms of two-parton GPDs. 

In the independent particle approximation which is used in all Monte Carlo models 
with multiparton interactions, the two-parton GPD is equal to the product of single 



particle GPDs discussed in section 2. Using parametrization of Eq. (10 1 one finds 

mm 

which leads to approximately a factor of two smaller cross section than the one observed 
at the Tevatron: TrR 2 nt 34m6 as compared to the experimental value of 7ri?? nt w 
15m6. The 3 — » 4 processes plays a minor role in the Tevatron kinematics and do not 
allow to solve this discrepancy. 

A fix implemented in the PYTHIA is to use a much more narrow distribution in p - 
effectively a dipole with m 2 — 2GeV 2 . This does decrease 7ri?? nt to the value observed 
at the Tevatron. However it is a factor of two smaller than the one determined from the 
analyses of the hard exclusive processes - see discussion in section 2 and in particular 
Fig. 3. Correspondingly in this model the difference between median p 2 for minimum- 
bias processes and processes with dijet trigger is a factor of 2 larger than what follows 
from the analysis of the HERA data, cf. Fig. [5j 

In principle there is a question of how well the separation of the 2 — > 4 processes 
was performed in the experimental studies. However the most recent DO analvsis|15| 
seems to indicate that the 2 —¥ 4 contribution in the kinematics used to determine 
7ri?f„ t is very small. 

This appears to leave us with only one possibility - presence of significant parton - 
parton correlations at a nonperturbative scale. Currently we are performing estimates 
of these correlations. We find that the correlations are indeed large and may explained 
the enhancement we discussed in this section. 
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5 Fluctuations of the gluon field and high multiplic- 
ity events at LHC 

Strength of the gluon field should depend on the size of the quark configurations. For 
example, the gluon field in the small configurations should be strongly screened - the 
gluon density much smaller than average. It is possible to extract from the comparison 
of the diffractive processes: 7^ + p — > V + X and 7^ + p — >• V + p the dispersion of the 
gluon strength at small x [15] : 



(G 2 ) - (G) 2 da T 



(G) 2 dt I dt 



(17) 

t=0 



The HERA data indicate that for Q 2 ~ 3GeV 2 and x ~ 10" 3 co g ~ 0.15 -7- 0.2 which is 
rather close to the value for the analogous ratio for the soft diffraction which measures 
fluctuations of overall strength of the soft hadronic interactions. 

How can one probe the gluon fluctuations in pp collisions? Let us consider multi- 
plicity of an inclusive hard process - dijet,... as a function of some cuts - for example 
overall hadron multiplicity: M (trigger) and build the ratio 

R = M (trigger) 

M (minimum bias) 

If there are no fluctuations of the parton densities, the maximal value of R is reached 
if the trigger selects collisions at small impact parameters 6^0. Using Eq. [I0]we find 

2 

777, 

R = P 2 (0)<7 in (pp) = -^o m (pp) » 4.5. (19) 

1Z7T 

Any larger enhancement of R could arise only from the fluctuations of the gluon density 
per unit area. 

The first measurement which maybe relevant for addressing the question of fluctu- 
ations was reported by ALICE [TS] . The multiplicity of Jj ip was studied as a function 
of the multiplicity in the central detector, namely 

dN ch d V = ,, M r , (20) 

{dNch/dri v= o} 

for dN^fJdr) < 5. It was found that R increases with increase of dN^/dr] reaching 
values fa 5 for dN^ h /dr) ~ 4. This number is close to what we estimated above. Any 
further increase of R would require presence of the fluctuations in transverse gluon 
density. An enhancement above the 6 = effect is given by the factor 

R ^ _ 9N(xi,Q 2 \n)g N (x 2 ,Q 2 \n) (S) ^ 
tl g N (xi 1 Q 2 )g N (x 2 ,Q 2 ) S 

Here n labels configurations selected by the trigger, and S is the area of the transverse 
overlap. In principle Rfi could reach very large values. For example, if we consider a 
collision of two protons in cigar shape configurations with the same gluon density for 
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different orientations of the protons, the enhancement would be proportional to the 
ratio of the principal axes of the ellipsoid. Another mechanism for the enhancement 



of Rfi is presence of the dispersion in the gluon density with uj g ~ 0.15 -j- 0.2 (Eq. 17) 
which leads to a few percent probability for the gluon field to be a factor 1.5 larger 
than average. 

These observations maybe of relevance for the discussion of the high multiplicity 
(HM) events studied by the CMS [T5]. In the analysis very rare events were selected 
which have the overall multiplicity for \rj\ < 2.4 of at least a factor of > 7 larger than 
the minimum-bias events. Probability of such events is very small: 

Phm « 10~ 5 - 10- 6 . (22) 

The two-particle correlations were measured as a function of the distance in the pseudo- 
rapidity - At] and the azimuthal angle - Ac/). Three types of correlations were observed : 
(a) very strong local correlation for Arj ~ 0, Ac/) ~ 0, (b) strong correlation for Ac/) ~ ir 
for a wide range of At?, (c) a weak correlation for 2 < \Ar]\ < 4.8, Ac/> ~ - so called 
ridge. 

The first question to address is how to get such a large multiplicity. It is pretty 
obvious that such events should originate from very central collisions. Based on our 
knowledge of P2Q)) we find that the probability of the collisions at b < 0.2 fm is ~2%. 
Using information about dispersion of fluctuations of the gluon fields we estimate the 
probability of fluctuation where both nucleons have g > 1.5giy(x) is > 10~ 3 . So a 
natural guess is that the CMS trigger selected central collisions with enhanced gluon 
fields in both nucleons. This should lead to a much higher rate of jet production per 
event. 

Indeed our analysis of the HM data indicates presence of a large total excess trans- 
verse momentum in the A</> ~ ir region: p° alance > 15GeV/c. Presumably it is due 
to production of two back to back jets with the trigger jet generating the narrow 
same side correlation. Qualitatively, a large probability of the dijets maybe due to the 
combination of centrality and the gluon density fluctuation. 

Note also that the increase of the multiplicity due to selection of b ~ and selection 
of b ~ and enhanced dijet production is not sufficient to generate a factor of 7 increase 
in the multiplicity - without of the gluon density fluctuations these two effects typically 
lead to N c h ~ 70. The g > gN( x ) gluon fluctuations would naturally lead to a further 
increase of N c h- 

The same side ridge could originate from the the string effect |12) . This could be 
tested by studying collisions with production of dijets with px ~ 15 GeV/c without 
HM trigger. Alternative mechanism would be fluctuations of the transverse shape of 
the colliding nucleons plus presence of the absorptive effects for p t < 3GeV/c. Such a 
scenario appears quite natural for the high density mechanism we discuss here. 



6 Where does the non-linear regime set in? 

In the leading log approximation one can derive a relation between the QCD evolution 
equations and the target rest frame picture of the interaction of small color dipoles 
with targets expressing it through the gluon density in the target, see |5D] and refs. 
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therein. Matching the behavior of the dipole cross section in the pQCD regime and of 
large size dipoles in the regime of soft interaction, it is possible to write interpolation 
formulae for the dipole - nucleon cross section for all dipole sizes and describe the total 
cross section of DIS at HERA. 

To determine how close is the interaction strength to the maximal allowed by the 
unitarity it is necessary to consider the amplitude of the qq dipole - nucleon interaction 
in the impact parameter space: 

r,?(*,&) = ^7^)2 / d 2 qe^A gS _ N (s,t), (23) 

where A qq -N(s, t) is the elastic amplitude of the qq dipole - nucleon scattering normalized 
to ImAqqpj = is q q-NO'tot(qq ~ N). The limit T qq -(s,b) = 1 corresponds to the regime 
of the complete absorption - black disk regime (BDR) - the maximal strength allowed 
by the S-channel unitarity. 

The t dependence of the qq dipole - nucleon elastic scattering amplitude can be 
obtained from the studies of the exclusive vector meson production in the regime where 
QCD factorization theorem for the exclusive processes allows to relate the t dependence 
of the amplitude to the t-dependence of the gluon GPDs. 

Combining this information with the information on the total cross section of the 
dipole - nucleon interaction allows us to determine T q q N (s, b) as function of the dipole 
size. A sample of the results for qq dipole -proton interaction which represent an 
update of the analysis of [U] is presented in Fig[7] For the case of color octet dipole 
Tfnei = (9/4)r^ ei leading to Y" much closer to one. As a result the gluon induced 
interactions are close to the BDR for a much larger range of the dipole sizes (this is 
consistent with the observation at HERA of a much larger probability of the diffraction 
in the gluon induce small x DIS processes). 

Note also that T = 1/2 already corresponds to a probability of inelastic interaction 
of 3/4 which is close to one. One can also demonstrate that the inelastic interactions 
get much larger corrections for the structure of the final states than the total cross 
section, see discussion in [2"2"] . 

In the case of the nuclear target the gain in the value of T(b ~ 0) is rather modest 
due to the leading twist shadowing. The main gain in the nucleus case is due to a weak 
dependence of T(6) on b for a broad range of b. 

Information about T q q-N(b) can be used to estimate the range of the transverse 
momenta for which probability of the inelastic interaction of parton is close to one [2] . 
The results of this analysis for 6 = are presented in FigjH] One can see from the figure 
that interaction of gluons is close to the BDR for a wide range of virtualities for the 
central pp collisions at the LHC. This is because a parton in the nucleon with a given 
x\ resolves the gluons in the second nucleon with x 2 down to Ap^/xis. For example 
taking x ~ 1(T 2 , </s=14 TeV and ~ p| = 4GeV 2 we find x 2 {min) = 1(T 4 . Though 
there are substantial uncertainties in this analysis due to the use of the leading log 
DGLAP approximation and extrapolation of the gluon densities to very small x, the 
analysis provides a reasonable estimate of the magnitude of pr (BDR) . Note also that 
the range we find for the RHIC kinematics for the central NA collisions is consistent 
with the effect of the suppression of the forward pion production due to the onset of 
the BDR which we discuss in Sects. 8, 9. 
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Figure 7: Impact parameter distribution of qq dipole interaction with protons adopted 
from pTj . 



Consequently the LT approximation should be broken for a wide range of x\, X2 for 
gluon - gluon interactions (these are x's for which DGLAP works well for the DIS at 
HERA). A breakdown of the LT approximation maybe of relevance for interpretation 
of the empirical observation made in a number of the Monte Carlo studies that to 
avoid contradictions with the data one needs to introduce a cutoff for minimal pt of 
the hard interactions. The cutoff is pT{min) > 3 GeV/c for y/s = 7 TeV and grows 
with s. In the Monte Carlo models the value of this cutoff is effectively driven by the 
requirement that the multiplicity of hadrons due to hard interactions remains below 
the total multiplicity observed experimentally. One can reach similar results (avoiding 
questions about sensitivity to the hadronization mechanism) based on the calculation of 
the probability of hard interactions as a function of the impact parameter and requiring 
that it does not exceed Ti ne i(b) which is known from the S-channel unitarity [23l 124) . 
Note however that the MC models with ppirnin) cutoff strongly reduce the interactions 
of the large x partons with the target which contradicts the proximity to the BDR. 
Hence one may expect that such Monte Carlo models overestimate the cross section of 
production of large xp hadrons, especially for the central collisions. 
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Figure 8: The p T range where interaction is close to the BDR for the interaction of qq 
and color octet dipoles plotted as a function of the energy of the dipole and of x of the 
interacting parton for pp interactions at y/s = 14TeV 



7 Postselection effect in BDR - effective fractional 
energy losses 

It was demonstrate in 25j that in the BDR interactions with the target select config- 
urations in the projectile wave function where the projectile's energy is split between 
constituents much more efficiently than in the DGLAP regime. The simplest example 
is inclusive production of the leading hadrons in DIS for Q < 2pr(BDB). Interac- 
tions with the target are not suppressed up to pr ~ pt(BDR), leading to selection of 
configurations in 7* where longitudinal fractions carried by quark and antiquark are 
comparable. The photon energy splits between the partons before the collision. It is 
the interaction that selects at different energies different set of configurations which are 
resolved. Hence we refer to this phenomenon as the postselection. As a result to a first 
approximation the leading hadrons are produced in the independent fragmentation of 
q and q: 

D^\z) = 2 ^ dyD h q {z/y)- A {\ + (2y - l) 2 ), (24) 
leading to a strong suppression of the hadron production at xp > 0.3. 
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In the case of a partem of a hadron projectile propagating through the nucleus near 
BDR effective energy losses were estimated in Ref.[5S]. For quarks they are expected 
to be of the order of 10 % in the regime of the onset of BDR and larger deep inside 
this regime. Also the effective energy losses are somewhat larger for gluons as the g — > 
gg splitting is more symmetric in the light cone fractions than the qg splitting. 

8 Leading hadron production in hadron - nucleus 
scattering 

Production of leading hadrons with p t ~ few GeV/c in hadron - nucleus scattering at 
high energies provides a sensitive test of the onset of the BDR dynamics. Indeed in this 
limit pQCD provides a good description the forward single inclusive pion production in 
pp scattering at RHIC [27 . At the same time it was found to overestimate grossly the 
cross section of the pion production in dAu collisions at RHIC in the same kinematics. 
The analysis of [25] has demonstrated that the dominant mechanism of the single pion 
production in the NN collisions in the kinematics studied at RHIC is scattering of 
leading quarks of the nucleon off the gluons of the target with the median value of x g 
for the gluons to be in the range x g ~ 0.01-^0.03 depending on the rapidity of the pion. 
The nuclear gluon density for such x is known to be close to the incoherent sum of 
the gluon fields of the individual nucleons since the coherent length in the interaction 
is rather modest for such x. As a result the leading twist nuclear shadowing effects 
can explain only a very small fraction of the observed suppression [25] and one needs a 
novel dynamical mechanism to suppress generation of pions in such collisions. It was 
pointed out in [28] that the energy fractional energy losses on the scale of 10% give a 
correct magnitude of suppression of the inclusive spectrum due to a steep fall of the 
cross section with xp which is consistent with the estimates within the postselection 
mechanism. 

An important additional information comes from the measurement of the corre- 
lation of the leading pion production with production of the pion production at the 
central rapidities [29j [30] . This corresponds to the kinematics which receives the dom- 
inant contribution from the scattering off gluons with x g ~ 0.01 -j- 0.02. The rate of 
the correlations for pp scattering is consistent with pQCD expectations. An extensive 
analysis performed in [26] has demonstrated that the strengths of "hard forward pion" 
- "hard r\ ~ pion" correlations in dAu and in pp scattering are similar. A rather small 
difference in the pedestal originates from the multiple soft collisions. Smallness of the 
increase of the soft pedestal as compared to pp collisions unambiguously demonstrates 
that the dominant source of the leading pions is the dAu scattering at large impact 
parameters. This conclusion is supported by the experimental observation [3T] that 
the associated multiplicity of soft hadrons in events with forward pion is a factor of 
two smaller than in the minimum-bias dAu events. A factor of two reduction factor 
is consistent with the estimate of [26] based on the analysis of the soft component of 
i] = production for the forward pion trigger. Overall these data indicate that (i) 
the dominant source of the forward pion production is the 2 — > 2 pQCD mechanism, 
(ii) production is dominated by projectile scattering at large impact parameters, (iii) 
proportion of small x g contribution in the inclusive rate is approximately the same for 
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pp and dAu collisions. 

A lack of additional suppression of the x g ~ 0.01 contribution to the double inclusive 
spectrum as compared to the suppression of the inclusive spectrum is explained in the 
post-selection mechanism as due to a relatively small momentum of the produced gluon 
in the nucleus rest frame putting it far away from the BDR. 

It is difficult to reconcile enumerated features of the forward pion production data 
with the 2 — ► 1 mechanism 52] inspired by the color glass condensate model. In the 
scenario of [32] incoherent 2 — > 2 mechanism is neglected, a strong suppression of the 
recoil pion production is predicted. Also it leads to a dominance of the central impact 
parameters and hence a larger multiplicity for the central hadron production in the 
events with the forward pion trigger. The observed experimental pattern indicates the 
models [55J which neglect contribution of the 2—^2 mechanism and consider only 2 — > 1 
processes strongly overestimates inclusive cross section due to the 2 — > 1 mechanism. 

Overall the observed regularities of inclusive forward pion production and forward 
central correlation phenomenon give a strong indication of breakdown of the pQCD 
factorization due to the propagation of high energy partons through the nuclear media. 
The modification of the nuclear gluon density at small x < 0.01 plays a small role in 
this kinematics. 

9 Production of two forward pions and double-parton 
mechanism in pp and dA scattering 

In Ref. [25] we suggested that in order to study the effects of small x gluon fields in 
the initial state one should study production of two leading pions in nucleon - nucleus 
collisions. Recently the data were taken on production of two forward pions in dAu 
collisions. The preliminary results of the studies of the reactions pp — > 7r°7r° + X, d — 
Au —} 7r°7T° + X, where one leading pion served as a trigger and the second leading 
pion had somewhat smaller longitudinal and transverse momenta [3U [55] . The data 
indicate a strong suppression of the back to back production of pions in the central 
dAu collisions. Also a large fraction of the double inclusive cross section is isotropic in 
the azimuthal angle Atp of the two pions. 

We performed a study of this process in [55] . We focused on explaining the isotropic 
component of the double pion spectrum and understand the origin of the suppression 
and in particular whether it is consistent with the post selection mechanism which we 
discussed above for the case of production of one forward pion. 

9.1 Forward dipion production in pp scattering 

It is instructive to start with the case of pp scattering. The leading twist contribution 
- 2 — > 2 mechanism - corresponds to the process in which a leading quark from the 
nucleon and a small x gluon from the target scatter to produce two jets with leading pi- 
ons. In this kinematics x g < M 2 (irir)/x q SNN- Production of two pions which together 
carry a large fraction of the nucleon momentum can occur only if x q is sufficiently close 
to one. The results of our calculation show that the average value of x q for typical cuts 
of the RHIC experiments is pretty close to one. 
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Obviously it is more likely that two rather than one quark in a nucleon carry 
together x close to one. This suggests that in the discussed RHIC kinematics pro- 
duction the "double-scattering" contribution with two separate hard interactions in a 
single pp collision could become important. Hence though the discussed contribution 
is a " higher- twist" , it is enhanced both by the probability of the relevant two quark 
configurations and the increase of the gluon density at small x which enters in the 
double-scattering in the second power. 

One can derive the expression for the double-scattering mechanism based on the 
analysis of the corresponding Feynman diagrams and express it through the double 
generalized parton densities in the nucleons which we discussed in section 4. Similar 
to Eq[l2] the cross section can be written in the form 



cfV 1 



dp T ,id m dp Ta d m irR mt abcda , b , c , d . 



/ dx a dxbdz c dx a 'dxb'dz c > 

J„/ „/ -J/ J 



^2 ^.ab^-cd *.a b —tc d 



^(...^(^(^W^W (25) 

Here f^, (x a , x a <) is the double parton distribution. If the partons are not correlated, 
it is equal to the product of the single parton distributions. For simplicity we neglected 
here correlations in the target as in our case x's for gluons are small. The dimensional 



factor nRf nt is given by Eq 15 In our numerical calculations we used ^R\ nt ~ 15 mb 
observed at the Tevatron which is smaller than the value obtained in the uncorrelated 
approximation (see discussion in section 4). 

We find that for the RHIC kinematics the only trivial correlation due to the fixed 
number of the valence quarks is important while the correlation between x a and x a > 
remains a small correction if we follow the quark counting rules to estimate the x a > 
dependence of f a ,a'(x a ,x a ') for fixed x a . The results of our calculation indicate that 
that the LT and double parton mechanisms are comparable for the kinematics of the 
RHIC experiments. This provides a natural explanation for the presence of a large 
component in the pp — > ir°ir° + X cross section measured in [311 which does not 
depend on the azimuth angle <p. In fact the number of events in the the pedestal is 
comparable to the number of events in the peak around <f> ~ 7r which is dominated by 
the LT contribution indicating that the LT and double-parton contributions are indeed 
comparable (see Fig. 10a below). 

Hence we conclude that the current experiments at RHIC have found a signal of 
double-parton interactions and that future experiments at RHIC will be able to obtain 
a unique information about double quark distributions in nucleons. It will be crucial 
for such studies to perform analyses for smaller bins of rj and preferably switch to the 
analysis in bins of Feynman x. 



9.2 Production of two forward pions and double-parton mech- 
anism in dAu scattering 

Let us extend now our results to the case of d-A scattering studied at RHIC. In this 
case there are three distinctive double-parton mechanisms depicted in Fig. [9j The 
first two are the same as in the pA scattering - scattering of two partons of the nucleon 
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Figure 9: Three double parton mechanisms of dipion production. 



off two partons belonging to different nucleons (mechanism a), and off two partons 
belonging to the same nucleon of the target (mechanism b) The third mechanism, 
which is not present for pA scattering is scattering of one parton of proton and one 
parton of the neutron off two partons of the nucleus. Let us consider the ratio of the 
double-parton and leading twist contributions for dA and pp collisions 

r dA = r a + r b + r c =°^r-^. (26) 
a LT {dA) oltXPV) 

The contribution to r^A of the mechanisms (a) , (c) is given by |38j : 

r c = T(b)a eff ;r a = l, (27) 

where T(b) is the standard nuclear profile function (J d 2 bT(b) — A). Here we neglected 
nuclear gluon shadowing effect which is a small correction for the double-parton mech- 
anism (cf. Ref.[2S]) but maybe important for the LT mechanism where x g maybe as 
low as 10 -3 due to the leading twist shadowing (see discussion below). For the cen- 
tral d — Au collisions Ta ~ 2.2 fm~ 2 and so r a /r c ~ 1/3. The contribution (b) can 
be calculated in a model independent way since no parton correlations enter in this 
case. The ratio of rj, and r c will be close to 1 at midrapidity, where correlations and 
valence-gluon scattering are not very important. Foward large rapidities, however, rj, 
must become much larger than r c , since it is not subject to the constraint x a + x a > < 1 
because of the fact that for (b) the proton and the neutron scatter independently. 

As a result TdAu for small b becomes of the order ten: r^A changes from ~ 9 to ~ 
12 for irR 2 int = 15 -7- 20 mb. 

Since the single inclusive pion spectrum for 772 ~ 2 -j- 3 is suppressed by a factor of 
the order RaQj) = 1/3 -j- 1/4 we find for the ratio of the pedestals in dAu and pp: 

Rpedestal = r A R A (b) ~ 2.5 -f- 4, (28) 

which should be compared with the experimental value of R pe destai ~ 3. Hence we 
naturally explain the magnitude of the enhancement of the pedestal in central dAu 



collision (see horizontal magenta lines in Fig 10 ) 

If most of the pedestal in the kinematics studied at RHIC is due to the double- 
parton mechanism, the uncertainties in the estimate of the rates due to this mechanism 
and uncertainties in the strength of the suppression of the single inclusive forward pion 
spectrum at b ~ would make it very difficult to subtract this contribution with a 
precision necessary to find out whether all pedestal is due to double-parton mechanism 
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Figure 10: STAR data for dipion production. Red curves are the Gaussian fit to the data. 
The horizontal magenta lines illustrate the strength of the double scattering while the solid 
black curve in the dA plot illustrates the effect of the reduction of the 2 — > 2 contribution 
by a factor of four as compared to the pp case. 

or there is a room for a small contribution of the 2 — > 1 broadening mechanism as 
it was assumed in [37]. Note also that in [37 authors calculated the ratio of the 
double inclusive cross section and the single inclusive cross section in the color glass 
condensate approach and compared this ratio with the data. However since the single 
inclusive spectrum is grossly overestimated by the model (see discussion in section 8) 
such procedure is not legitimate. 

The suppression of the away peak originating from the LT contribution is due to 
two effects: (i) the gluon shadowing for x ~ 10~ 3 and b < 3fm and Q 2 ~ few GeV 2 
reduces the cross section by a factor of about two, (ii) stronger effect of effective 
fractional energy losses due to larger x of the quark in the LT mechanism than in the 
double parton mechanism, leading to a suppression factor of the order two |36j . All 
together this gives a suppression of the order of four as compared to the single pion 
trigger, which is consistent with the STAR observation - see solid curve in Fig. [Tup . 
It corresponds to the overall suppression of the order of ten. This is pretty close to 
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the low limit for the suppression estimated as the probability of the "punch through" 
mechanism - contribution from the process where a quark scatters off one nucleon but 
does not encounter any extra nucleons at its impact parameter. Probability of such 
collisions at b ~ for interaction with Au nucleus is of the order 5 -j- 10% [39] . 

The data are consistent with suppression of the away peak by a factor > 4 and 
the reduction of the away peak relative to pedestal of the order of ten. The data 
may indicate that in addition to overall suppression there is some broadening of the 
away peak. Such effect is present in the postselection mechanism, though for the very 
forward kinematics it is a correction to the effective energy losses. 

For the LHC kinematics the discussed effects will be grossly amplified and extend 
to much wider range of x - the same parton - target energy corresponds to rescaling 
of x of the factor of slhc / s rhic > 10 3 - In addition for x < 0.1 the gluons give the 
dominant contribution while the BDR scale of p\ is about a factor of two larger in this 
case (cf. Figj|]). 



10 Conclusions 

Studies aimed at understanding the underlying dynamics of pp scattering at the LHC 
energies face a number of challenges. The challenges discussed in this lecture include 

1. Building models of inelastic collisions with realistic transverse parton distribu- 
tions. 

2. Including effects of the correlations between the partons in order to describe the 
rate of multiparton interactions. 

3. Realistic modeling of the BDR effects at moderate transverse momenta 

4. Describing forward production at the LHC which is most sensitive to the BDR 
dynamics and in particular to the effect of effective fractional energy losses. 
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